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Vcalc:  a  3-Space  Vector  Calculator 


INTRODUCTION 

Vcalc  is  a  3-space  vector  calculator  providing  graphic  displays  for  vector  and  scalar  operations. 
Designed  to  operate  on  a  Sun  Workstation,  the  program  may  be  run  stand-alone  or  as  a  debugging 
tool  linked  to  a  C  program.  It  is  a  straightforward  and  useful  tool  whose  development  was  moti¬ 
vated  by  the  difficulty  of  debugging  geometric  programs. 

STAND-ALONE  MODE 

To  run  Vcalc  stand-alone,  simply  type  vcalc.  The  screen  as  shown  in  Figure  1  will  appear. 


Vcalc  COMMANDS 


To  select  an  operation,  type  the  corresponding  acronym.  For  example,  to  receive  instruc¬ 
tions  for  computing  the  scalar  triple  product  of  three  vectors,  type  tp.  Following  is  a  brief  descrip¬ 
tion  of  each  command: 


as  Store  the- sum  of  two  scalars  into  the  output  scalar  slot, 

av  Store  the  sum  of  two  vectors  into  the  output  vector  slot, 

cl  Qear  all  vector  and  scalar  slots. 

cp  Store  the  cross  prcxiuct  of  two  vectors  into  the  output  vector  slot. 

cs  Copy  a  new  scalar  into  a  scalar  slot,  or  copy  an  existing  scalar  from  one  slot  into 

another. 

cv  Copy  a  new  vector  into  a  vector  slot,  or  copy  an  existing  vector  from  one  slot  into 

another. 


d 

dp 

ds 

dsqr 


Store  the  distance  between  two  vectors  into  the  output  scalar  slot. 

Store  the  dot  product  of  two  vectors  into  the  output  scalar  slot. 

Store  the  quotient  of  two  scalars  into  the  output  scalar  slot 

Store  the  distance  squared  between  two  vectors  into  the  output  scalar  slot. 
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ins  Store  the  product  of  two  scalars  into  the  output  scalar  slot 

mvs  Store  the  product  of  a  vector  and  scalar  into  the  output  vector  slot, 

nv  Store  a  normalized  vector  into  a  specified  output  vector  slot, 

q  Exit  Vcalc 

rv  Rotate  a  vector  N  degrees  around  another  vector  and  store  the  result  into  a  specified  output 

slot. 

rx  Rotate  a  vector  N  degrees  around  the  X  axis  and  store  the  result  into  a  specified  output 

slot. 

ry  Rotate  a  vector  N  degrees  around  the  Y  axis  and  store  the  result  into  a  specified  output 

slot. 

rz  Rotate  a  vector  N  degrees  around  the  Z  axis  and  store  the  result  into  a  specified  output  slot. 

sc  Adjust  the  graduation  spacing  (i.e  the  scale  of  the  graph). 

ss  Store  the  difference  between  two  scalars  into  the  output  scalar  slot. 

sv  Store  the  difference  between  two  vectors  into  the  output  vector  slot. 

tp  Store  the  scalar  triple  product  of  three  vectors  into  the  output  scalar  slot. 

vm  Compute  the  length  of  a  vector  and  store  the  result  into  the  output  scalar  slot. 

vsq  Store  the  dot  product  of  a  vector  with  itself  into  the  output  scalar  slot. 

X  View  the  vector  display  along  the  X  axis. 

y  View  the  vector  display  along  the  Y  axis. 

z  View  the  vector  display  along  the  Z  axis. 

Whereas  several  vector  operations  allow  the  user  to  specify  in  which  “Vector  List”  slot  the 
result  is  to  be  placed,  the  output  vector  will  always  appear  in  slot  0  (marked  “Output  Vector”)  by 
default.  T  ame  concept  applies  to  scalars. 

DEBUGGING  C  PROGRAMS 

To  run  Vcalc  as  a  debugging  tool,  two  source  files  need  to  be  included  within  the  compilation 
stream,  vcalc.c  and  vcjib.c.  Pixrects  and  Math  are  the  only  libraries  required.  An  example  compila¬ 
tion  might  look  like; 
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cc  myprogx  vcalcc  vcjib.c  >0  myprog  'Im  'Ipixrect 

Vcalc  involves  the  use  of  the  two  macros  ev(vector),  es(scalar),  and  the  function  vc(): 
finclude  <vcalc.h> 

ev(vcctor)  /*  ev  is  used  to  load  a  vector  into  Vcalc.*/ 

double  *vector,  /*  The  input  argument  points  to  an  array  of  three  doubles.  */ 

es(scalar)  /*  es  is  used  to  load  a  scalar  into  Vcalc.*/ 

double  scalar,  /*  The  input  argument  contains  a  scalar.  */ 

vc();  /*  vc()  activates  Vcalc  along  with  all  vectors  and  scalars 

previously  loaded  with  ev()  and  es().  */ 

The  following  program  provides  a  simple  example  of  how  Vcalc  is  called  from  a  C  program. 
Note  the  use  of  the  above  macros  in  poning  program  data  into  Vcalc. 

#include  <stdio.h> 

#include  “vcalc.h” 

double  vecl[3],  vec2[31,  vec3[31.vec_array|  1 1113]; 

double  scall,  scal2[12]; 

mainO 

(  ■  - 

veclI0]=3.3;  vecl|  1 1=4.4;  vccl(21=5.5; 

ev(vecl);  /‘Add  vecl  to  the  (currently  empty)  vector  list  */ 

vec2I0]=23.273;  vec2|  1)=-10.4;  vcc2121=15; 

ev(vec2);  /•  Add  vec2  to  the  vector  list  */ 

vec_arrayI5][0]=40; 

vec_aiTay[5][  1 1=9.K9, 

vec_array[5][2)=-20.4 1 ; 

cv(vec_array(5));  /•  Add  vcc_anay|5]  to  the  vector  list  */ 

vec3[0]=-28.00;  vcc.3|  1|=  3K;  vec3I2)=4.00; 

ev(vec3);  /•  Add  vec3  to  the  vector  list  */ 

seal  1=1 23.987; 

scal2[6]=3; 

es(scall);  /•  Add  scall  to  the  (currently  empty)  scalar  list  */ 

es(scal2[6]);  /•  Add  scal2  to  the  scalar  list  */ 

vc();  /•  Bring  up  Vcalc  */ 

) 

This  example  results  in  the  screen  image  shown  in  Figure  2.  The  figure  also  illustrates  the 
selection  of  the  rz  command.  Recall  that  the  rz  command  rotates  a  vector  around  the  Z  axis. 
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The  four  vectors  and  the  two  scalars  passed  to  vc()  via  ev()  and  es()  are  listed  in  the  “Vectoi 
List”  and  “Scalar  list.”  respectively.  The  vectors  are  drawn  and  labeled  according  to  their  index  into 
the  vector  list.  The  standard  vector  input  prompt  is: 

enter  vector  (0  thru  10  -or- 11  for  a  new  vector): 

0-10  refers  to  the  vector  list  index.  A  new  user-specified  vector  may  be  entered  by  selecting 
11.  Vector  number  4  has  been  selected  for  this  rotation.  In  our  test  program,  the  variable  name  vec3  is 
displayed  in  vector  slot  4  along  with  its  X,Y,Z  values.  Vector  slot  6  has  been  selected  as  the  output 
location  and  a  90  degree  rotation  is  specified.  The  result  of  rotating  (-28,-38,4)  90  degrees  around  the 
Z  axis  is  shown  in  Figure  3. 

Notice  in  Figure  3  that  the  output  from  the  previous  rotation  operation  was  stored  in  vector 
slot  6  without  a  name.  Only  those  vectors  and  scalars  passed  directly  from  the  calling  program  have 
their  names  listed  (vecl,  scal2  [6],  etc.).  Output  values  from  Vcalc  operations  as  well  as  vectors  and 
scalars  entered  in  by  hand  are  listed  without  variable  names. 

Figure  3  shows  that  we  have  selected  the  rotate  vector  command  (rv)  for  our  next  example. 
Vector  2  will  be  rotated  around  vector  1  by  180  degrees  and  the  output  will  be  stored  in  vector  2.  Vari¬ 
able  vec2  will  hold  the  result  of  rotating  the  previous  contents  of  vec2  180  degrees  around  vecl. 

The  current  graduation  spacing,  as  shown  directly  beneath  the  graph,  indicates  that  each  line 
on  the  graph  represents  one  unit  of  measure.  This  value  is  adjustable.  For  example,  we  can  set  each 
line  equal  to  five  units  by  simply  selecting  adjust  scaling  (as)  and  entering  5  at  the  prompt  as  shown 
in  Figure  4.The  resulting  screen  is  shown  in  Figure  5.  Figure  5  also  shows  us  instructing  Vcalc  to  nor¬ 
malize  vector  4  (i.e.  vec3)  and  place  the  result  in  vector  slot  8. 

In  Figure  6  we  take  the  cross  product  of  vector  8  (i.e.,  the  output  of  the  previous  operation) 
and  vector  I.  The  result  is  stored  in  the  default  Output  Vector  slot  (slot  0).  Besides  the  cross  product. 
Figure  7  shows  that  we  have  set  the  graduation  spacing  back  to  1.  We  can  view  the  current  vectors 
along  the  Y  axis  by  entering  Y  (Fig.  8). 

Several  scalar  operations  are  available.  Figures  9  and  10  show  us  taking  the  dot  product  of 
vectors  3  and  1,  with  the  output  going  into  the  default  output  scalar  slot  (scalar  slot  0). 

As  with  vectors,  we  always  have  the  option  of  using  new  scalar  data  by  entering  1 1  at  the  fol¬ 
lowing  prompt: 

enter  scalar  (0  thru  10  -or- 11  for  a  new  scalar): 

We  will  expand  the  capabilities  of  Vcalc  as  the  need  arises.  Since  one  of  our  major  interests  is 
the  interpretation  of  range  images,  one  possible  direction  for  expansion  is  the  fitting  of  simple  sur¬ 
faces  to  range  data,  and  the  computation  of  angular  and  linear  displacements  between  such  structures 
as  planes,  lines,  dihedral  edges,  trihedral  vertices,  etc. 
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Figure  1 .  Vcalc’s  opening  screen  -  In  this  example,  no  vector  or  scalar  data  has  been  passed  to  Vcalc. 


Figure  2.  The  user  is  instructing  Vcalc  to  rotate  vec3  90  degrees  around  the  Z  axis,  and  place  the  result  in  vector 
slot  number  6. 


Figure  3.  The  user  is  instructing  Vcalc  to  rotate  vec2  180  degrees  around  vecl,  and  place  the  result  in  vector  slot 
number  2. 
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Figure  4.  Vector  number  2  (vec2)  has  been  rotated  180  degrees  around  vecl.  The  graduation  spacing  (scaling)  is 
being  adjusted  from  1.0  to  5.0. 


•sic 


Figure  5.  With  graduation  spacing  now  equal  to  5.0,  Vcalc  is  being  instructed  to  normalize  vector  number  4  (vec3). 
The  output  will  be  placed  in  vector  slot  8. 
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The  user  is  instructing  Vcalc  to  take  the  cross  product  of  vector  number  8  and  vector 
The  output  will  be  automatically  stored  in  vector  slot  0. 
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Figure  7.  Graduation  spacing  has  been  set  back  to  1.0,  and  the  cross  product  of  vector  number  8  and  vector  number 
one  (vecl)  is  now  stored  in  vector  slot  0  as  well  as  being  graphically  displayed. 
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Figure  8.  By  selecting  Y,  we  can  view  the  cunent  vector  list  along  the  Y  axis. 
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Figure  9.  By  having  selected  Z  (not  shown),  we  are  once  again  viewing  the  vector  list  along  the  Z  axis.  The  user  is 
now  instructing  Vcalc  to  compute  the  dot  product  of  vector  number  3  (vec_array[S])  and  vector  number  1 
(vecl).  The  result  will  be  stoi^  in  scalar  slot  0. 
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Figure  10.  The  dot  product  of  vec_array(51  and  vcci  is  now  stored  in  scalar  slot  0. 


